Steady state thermodynamics for heat conduction 
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Following the proposal of steady state thermodynamics (SST) by Oono and Paniconi, we develop a 
phenomenological theory for steady nonequilibrium states in systems with heat conduction. We find 
that there is essentially a unique consistent thermodynamics, and make concrete predictions, i.e, the 
existence of a new osmotic pressure and a shift in the coexistence temperature. These predictions 
allow one to test for the quantitative validity of SST by comparing them with experiments. 
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Construction of a statistical mechanics that apply to 
nonequilibrium states has been a challenging open prob- 
lem in theoretical physics But so far it is not known 
how the desired probability measures for nonequilibrium 
states look like, or even whether the measures can be 
written in compact forms |^ . Recalling the history that 
the conventional thermodynamics was an essential guide 
when Boltzmann, Gibbs, and others constructed equilib- 
rium statistical mechanics, it may be a good idea to start 
from the level of thermodynamics. 

The standard theory of nonequilibrium thermodynam- 
ics is based on local equilibrium hypothesis, which 
roughly asserts that each small part of a nonequilibrium 
state can be regarded as a copy of a suitable equilib- 
rium state. But such a description seems insufficient for 
general nonequilibrium states. Consider, for example, a 
system with a steady heat flow. It is true that the quan- 
tities like the temperature and the density become essen- 
tially constant within a sufficiently small portion of the 
system. But no matter how small the portion is, there 
always exists a heat flux passing through it and hence the 
local state is not isotropic. This suggests that the local 
state cannot be identical to an equilibrium state (which 
is isotropic), but should be described rather as a local 
steady state. 

Among existing attempts in nonequilibrium thermody- 
namics to go beyond local equilibrium treatments , the 
steady state thermodynamics (SST) proposed by Oono 
and Paniconi Q seems to be most sophisticated and 
promising. The basic strategy of , in our own interpre- 
tation, is to i) look for a thermodynamics which describes 
a steady state as a whole, ii) clarify operational pro- 
cedures for determining thermodynamic quantities, and 
iii) respect the general mathematical structure of thermo- 
dynamics. In the present Letter, we apply this strategy 
to a concrete problem of heat conduction (in a fluid), 
and show that theoretical consistency leads one to an es- 
sentially unique thermodynamics. We then make some 
concrete predictions which may be confirmed quantita- 
tively in experiments. Extensions to other systems and 
discussions of related microscopic results will appear in 
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FIG. 1: A typical system treated in SST. The left and the 
right walls have temperatures T and T' , respectively, and 
there is a steady heat flux J. We get a local steady state by 
restricting our attention to the thin portion within a distance 
L from the left wall (as denoted by the dashed line) 



Local steady state: We consider a macroscopic system 
of a single substance in a closed cylindrical container with 
the cross section area A. The left and right walls of the 
container have very efficient heat conduction, and are 
kept at constant temperatures T and T", respectively, 
with the aid of external heat baths. See Fig. |l|. The side 
walls of the container are perfectly adiabatic. 

If the system is kept in this setting for a sufflciently 
long time, it is expected to reach a unique steady state 
with a constant heat current without any macroscopically 
observable changes. (We assume that convection does 
not take place.) By J we denote the total energy that 
flow into the system from the right wall within a unit 
time. 

We now restrict our attention to a thin portion of the 
system within a small distance L from the left wall as in 
Fig. |l]. Here the length L is taken so that the tempera- 
ture in the thin portion becomes essentially constant. We 
assume (as in local equilibrium approaches) this is real- 
ized with L which is much larger than any microscopic 
scales. The state in this thin portion is the local steady 
state that we study 0|. Our first crucial assumption is 
that the local steady state can be fully specified by four 
macroscopic parameters as (T, J; V, N) , where V = AL 
is the volume and A'' is the amount of substance in the 
thin portion. We stress that the heat flux J need not be 
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FIG. 2: The porous wall in the middle of the container sepa- 
rates an equilibrium state and a steady nonequilibrium state. 
If we restrict our attention to the vicinity of the porous wall 
(denoted by dashed lines), we get a situation where a local 
steady state and an equilibrium state are in balance with each 
other. We use this setting to measure the chemical potential 
of the steady state. We will also see that there appears a 
force (of nonequilibrium origin) that pushes the porous wall 
towards the equilibrium region irrespectively to the sign of 
J 7^ 0. This is an example of "flux induced osmosis." 



small. 

As in the conventional thermodynamics, it is essential 
to consider decomposition/recombination and scaling of 
local steady states. In doing so, we shall always fix the 
cross section A and vary only the length L (within the 
range the system remains thin). Consider splitting the 
system of length L into those with lengths Li and L2 
(with L = Li + L2) hy e. plane parallel to the left and 
right walls Q . Correspondingly, we assume that the local 
steady state (T, J; V, N) can be decomposed into two lo- 
cal steady states (T, J; Vi, iVi) and (T, J; V2, A^2), and the 
two states can be reconibined back into (T, J; V, N). Here 
Vi = ALi, and Ni+N2 = N. Similarly, for A > 0, we as- 
sume that one can scale the length L to XL to get a scaled 
copy (T, J; XV, XN) of the state (T, J; V, N). These ob- 
servations imply that the heat flux J behaves in a similar 
way as T, and hence should be regarded as an "intensive" 
variable. This identification is essential in our theory. 

SST free energy: Our second essential assumption is 
the existence of the SST free energy F{T, J; V, N), which 
is a thermodynamic potential describing the response of 
the local steady state when the extensive variables V 
and N are varied. We assume that there is a func- 
tion F(T, J; V, N) which is concave in intensive vari- 
ables T, J, and convex in extensive variables V , N . It 
should have additivity F{T, J; V, N) = F{T, J; Vi, iVi) + 
F{T, J;V2, N2) for any decomposition as above, exten- 
sivity F{T,J;XV,XN) = XF{T,J;V,N) for A > 0, and 
symmetry F{T, J; V, N) = F{T, - J; V, N) to reflect the 
obvious left-right symmetry of thin systems. Most impor- 
tantly, we require that the derivatives of F{T, J; V, N) in 
V and N have operational meanings exactly parallel to 
those in the conventional thermodynamics. 

More precisely, we first require dF{T, J;V, N)/dV = 
—p{T, J; V, N) , where the pressure p{T, J; V, N) is sim- 
ply determined by measuring the force that the system 
exerts on the left (or the right) wall 0. We then require 



dF{T, J; V, N)/dN = n{T, J; V, N), where the existence 
and measurability of the chemical potential /x(T, J; V, N) 
are assumed. A device for measuring fj,(T, J; V, N) 
is depicted in Fig. ||. Here the left and the right walls of 
the container are kept at constant temperatures T and 
T' , respectively. In the middle of the container, there is 
another wall made of a porous medium which is kept at a 
constant temperature T with the aid of an external heat 
bath. We assume that the substance can move across the 
porous wall. Suppose that the whole system has reached 
a steady state. Then the left half of the system has a con- 
stant temperature, and is in a normal equilibrium state 
(provided that the middle walls is in a very efficient con- 
tact with a bath). The right half is in a nonequilibrium 
steady state (which is not necessarily local) with a con- 
stant heat flux. In order to examine the balance between 
the two parts, we restrict our attention to a thin part of 
the system within a fixed distance from the porous wall as 
denoted by dashed lines in Fig. |2[ If the distance is small 
enough, the part of the steady state can be regarded as 
local. Then we get a situation where a local steady state 
(T, J; V, N) and an equifibrium state (T, 0; V , N') are in 
balance with each other with respect to the exchange of 
the substance. Then it is natural to define /x(T, J; V, N) 
to be equal to ^{T,0;V' , N'), where the latter can be 
determined within the conventional equilibrium thermo- 
dynamics. The two derivatives and the extensivity de- 
termine F{T, J; y, N) without any ambiguities |0, as is 
obvious from the Euler equation 

F(T, J; V, N) = -Vp{T, J; V, N) + N ^l{T, J; V, N), 

which is derived by using the extensivity as usual. 
Let us define (extensive) SST entropy as 

dF{T,J;V,N) 



(1) 



S{T, J;V,N) 



dT 



and a new "extensive" quantity 

9F(T, J; V, N) 



^{T,J;V,N) 



dJ 



(2) 



(3) 



which we shall call nonequilibrium order parameter. 
From the symmetry and the concavity of F, one finds 
that ^{T,J;V,N) = -*(r, - J; A^) > for J > 0. 
We expect to have ^'(r, J; A^) > for J > in generic 
systems. Note that both S and 4" can be measured op- 
erationally since F can be. As in the conventional ther- 
modynamics, one can derive various identities between 
thermodynamic quantities. 

The existence of the SST free energy with the desired 
properties is nothing more than an optimistic assump- 
tion. It is possible in principle that such a theoretical 
framework as SST simply does not exist in Nature. We 
therefore predict two concrete phenomena — the exis- 
tence of a new osmotic pressure called flux-induced os- 
mosis (FIO) and a shift of coexistence temperature — 
and present some exact relations, which enable one to 
test for the quantitative validity of SST through experi- 
ments. 
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Flux-induced osmosis: Let us consider the situation 
in Fig. ^ and examine the behavior of the pressure 
Pss — p{T, J; V, N) of the steady state. We fix the tem- 
perature T and the pressure Pcq = p{T,0;V' , N') (by 
suitably varying the volume V') of the equilibrium state, 
and vary only the flux J. The chemical potential of the 
equilibrium state is thus constant, and so is the chemi- 
cal potential fi{T, J; V, N) of the local steady state (by 
definition). We differentiate the Euler equation (|]) by J 
with T and Peq fixed. Keeping in mind that V and N 
may depend on J, we find 



9p,s(r,Peq,J) _ ^{T,J;V,N) 



dJ 



V 



(4) 



Noting that Pss{T,Peq, 0) = Pcq and recalling the sign of 
^, this implies Pss > Poq in general. If ^' is nonvanishing 
(as we expect) then one has Pss > Pcq for J ^ 0. Recall- 
ing that the pressures are defined from mechanical forces 
exerted on the walls, this implies that the porous wall is 
actually pushed towards the equilibrium region irrespec- 
tively to the sign of the heat flux J. We stress that this 
force, which is absent in equilibrium or local-e quil ibrium 
treatments, is of purely nonequilibrium origin 12 . This 
is an example of a general phenomenon that we call flux- 
induced osmosis (FIO) ||]. 

In the same situation, one can also derive Q a non- 
trivial identity 



dpss{T,Pcq, J) _ Veq 
dpcq Vss ' 



(5) 



where v^q = V /N' and Vss = V/N. Since the identity 
involves only directly measurable quantities, it may be 
useful in quantitative tests of SST in (real and numerical) 
experiments. 

Shift of coexistence temperature: Consider again the 
original setting in Fig. |l|, and assume that the pressure 
p is kept constant. We further assume that a phase co- 
existence takes place in the system, i.e., the lower tem- 
perature region of the container is occupied by one phase 
(e.g. liquid) while the higher temperature region by an- 
other (e.g. gas). We then ask what is the temperature 
Tc{p, J) at the boundary between the two phases. Within 
local equilibrium treatments, one simply concludes that 
Tc(p, J) is the same as its equilibrium value Tc{p,0). In 
SST, however, we find (from an analysis similar to that 
for FIO) that 



dTcjp, J) _ V'high - '0: 



low 



dJ 



^high -^low 



(6) 



were V'low and "iAhigh are the molar ^! of the two phases 
at the coexistence point, and siow and Shigh are the cor- 
responding molar SST entropies. The identity (^) means 
that in general Tc(p, J) shifts from its equilibrium value. 
In case Tc(p, J) > Tc{p,0) one has a remarkable 

phenomenon of "heat flux induced freezing", i.e., one ob- 
serves a solid phase in a system with one wall having a 
temperature slightly higher than the melting point and 
the other wall having a much higher temperature. 



Choice of nonequilibrium variable: Finally let us 
make an important remark about the choice of nonequi- 
librium thermodynamic variable. A local steady state 
may also be specified as {T;V,N,t), where r is the 
(small) temperature difference between the left and right 
ends of the system. Since r is proportional to the length 
L in a thin uniform state, it should be regarded as 
an "extensive" quantity. Then, in any mathematically 
"healthy" thermodynamics, the corresponding free en- 
ergy F{T; V,N,t) (if exists) should be convex in V, N, 
and r. 

Following the standard argument, the convexity and 
the extensivity implies a variational principle 

F{T; Vi,N,,T*)+ F{T; F2, iVa, r - r*) 

= min{F(T; V^i, iVi, r') + F(T; V2,N2, r - r')}, (7) 

t' 

for fixed Vi, V2, Ni, N2, and r. This relation de- 
termines the temperature T -\- t* at the boundary of 
the two parts of the system (with volumes Vi and 
V2, respectively) when the total temperature difference 
T is fixed. In a differential form, this condition be- 
comes v{T;Vi,Ni,T*) = v{T;V2,N2,T - T*), where 
i^{T; V, N, t) = dF{T; V, N, T)ldT. 

On the other hand, we already know from the energy 
conservation law that r* can be determined by the con- 
dition J(r; Vi, iVi, r*) = J(r; F2, iVa, r - r*), where we 
expressed the heat flux J as a function of {T\V,N ,t). 
Since the conditions written in v and in J must be equiv- 
alent, we conclude that there is a function / such that 
J(T; V, N, t) = f{v{T- V, N, r)) for any (T; V, N, r). But, 
noting that J has a dimension of energy divided by time 
while V is dimensionless, we find that such a universal 
function / simply does not exist. This observation im- 
plies that one can never get a consistent thermodynamics 
by using the temperature difference t as a nonequilibrium 
variable. In fact we believe that the representation in 
terms of (T, J; V, N), where the nonequilibrium variable 
J directly reflects the energy conservation law, provides 
an essentially unique consistent thermodynamics for a 
system with steady heat flux [Q. 

Discussions: By following the philosophy of we 
are led to an essentially unique thermodynamics (SST) 
for heat conduction. The uniqueness suggests that, if a 
consistent thermodynamics for steady heat conduction 
exists at all, then it should be equivalent to what we 
have described here. The resulting theory led us to novel 
predictions which allow one to test for the quantitative 
validity of SST through experiments, other theoretical 
results 
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Let us stress that SST is by no means incompatible 
with the long range (spatial) correlation gencrically found 
in nonequilibrium steady states |l7|| . The long range cor- 
relation comes from the global description based on fluc- 
tuating hydrodynamics, while local aspects are conven- 
tionally described as local steady states. The aim of SST 
is to provide a more sophisticated alternative of the lat- 
ter, i.e., to give a unified description of nonequilibrium 
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modifications of the (local) equations of states. Since 
there is a big separation in the length scales, the exis- 
tence of the long range correlation does not conflict with 
the validity of SST. 

Perhaps the most crucial point about the potential sig- 
nificance of SST is whether it becomes a useful guide 
in the (future) construction of statistical mechanics for 
steady nonequilibrium states. The fact that we have ar- 
rived at an essentially unique theory is rather encourag- 
ing. We hope that, by trying to construct a statistical 



theory that is consistent with the (unique) nonequilib- 
rium thermodynamics, we are naturally led to a meaning- 
ful and correct statistical mechanics for steady nonequi- 
librium states. 
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